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Abstract 

The first observation of the decay — t- D^K*^ using pp data collected by the 
LHCb detector at a centre-of-mass energy of 7 TeV, corresponding to an integrated 
luminosity of 36 pb~^, is reported. A signal of 34.4 =b 6.8 events is obtained and the 
absence of signal is rejected with a statistical significance of more than nine standard 
deviations. The — )■ D^K*^ branching fraction is measured relative to that of 

W D^p^: ^g||r^o^o)^ = 1-48 ± 0.34 ± 0.15 ± 0.12, where the first uncertainty 
is statistical, the second systematic and the third is due to the uncertainty on the 



ratio of the and B^ hadronisation fractions. 
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1 Introduction 



A theoretically clean extraction of the Cabibbo-Kobayashi-Maskawa (CKM) unitarity 
triangle angle 7 can be performed using time-integrated B — )■ DX decays by exploiting 
the interference between Cabibbo-suppressed 6— )■ u and Cabibbo-allowed 6—)- c transitions 
[1, 2, 3, 4, 5, 6]. One of the most promising channels for this purpose is — )■ DK*^, 
where D represents a or a meson. ^ Although this channel involves the decay of a 
neutral B meson, the final state is self-tagged by the flavour of the K*^ so that a time- 
dependent analysis is not required. In the B^ — )■ DK*^ decay, both the B'^ — )■ D^K*^ 
and the B^ — )■ D^K*^ are colour suppressed. Therefore, although the B^ DK*^ decay 
has a lower branching fraction compared to the B^ DK^ mode, it could exhibits an 
enhanced interference. 

The Cabibbo-allowed S° D^K*^ and 5° d*^K*^ decays potentially provide a 
significant background to the Cabibbo-suppressed — )■ D^K*^ decay. The expected size 
of this background is unknown, since the -B°— )■ D^*^^K*^ decay has not yet been observed. 
In addition, a measurement of the branching fraction of B^ — D^K*^ is of interest as a 
probe of SU(3) breaking in colour suppressed B^^ — )■ D^V decays [7, 8], where V denotes 

a neutral vector meson. Thus, the detailed study of B^ — t- D^K*^ is an important goal 
with the first LHCb data. 

The LHCb detector [9] is a forward spectrometer constructed to measure decays of 
hadrons containing b and c quarks. The detector elements, placed along the collision axis 
of the Large Hadron Collider (LHC), start with the Vertex Locator, a silicon strip de- 
vice that surrounds the pp interaction region with its innermost sensitive part positioned 
8 mm from the beam. It precisely determines the locations of the primary pp interac- 
tion vertices, the locations of the decay vertices of long-lived hadrons, and contributes 
to the measurement of track momenta. Other tracking detectors include a large-area 
silicon strip detector located upstream of the 4 Tm dipole magnet and a combination of 
silicon strip detectors and straw drift chambers placed downstream. Two Ring-Imaging 
Cherenkov (RICH) detectors are used to identify charged hadrons. Further downstream 
an electromagnetic calorimeter is used for photon detection and electron identification, 
followed by a hadron calorimeter and a muon system consisting of alternating layers of 
iron and gaseous chambers. LHCb operates a two stage trigger system. In the first stage 
hardware trigger the rate is reduced from the visible interaction rate to about 1 MHz 
using information from the calorimeters and muon system. In the second stage software 
trigger the rate is further reduced to 2 kHz by performing a set of channel specific selec- 
tions based upon a full event reconstruction. During the 2010 data taking period, several 
trigger configurations were used for both stages in order to cope with the varying beam 
conditions. 

The results reported here uses 36 pb~^ of pp data collected at the LHC at a centre- 
of-mass energy ^/s = 7 TeV in 2010. The strategy of the analysis is to measure a ratio 
of branching fractions in which most of the potentially large systematic uncertainties 

^In this Letter the mention of a decay will refer also to its eharge-conjugatc state. 
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cancel. The decay — )■ D^p^ is used as the normahsation channeL In both decay 
channels, the is reconstructed in the Cabibbo-allowed decay mode ^ K~tt~^; the 
contribution from the doubly Cabibbo-suppressed — K'^ir^ decay is negligible. The 
K*'^ is reconstructed in the K*'^ K~^7r~ decay mode and the p° in the p°— vt+tt" decay 
mode. The main systematic uncertainties arise from the different particle identification 
requirements and the pollution of the 5° — ?• D'^p^ peak by 5° — )■ D°7r+7r~ decays where 
the 7r~^7T~ pairs do not originate from a p° resonance. In addition, the normalisation of the 
decay to a B^ decay suffers from a systematic uncertainty of 8 % due to the current 
knowledge of the ratio of the fragmentation fractions fs/ fd = 0.2671qo2o [10]- 

2 Events selection 

Monte Carlo samples of signal and background events are used to optimize the signal 
selection and to parametrize the probability density functions (PDFs) used in the fit. 
Proton beam collisions are generated with PYTHIA [11] and decays of hadronic particles 
are provided by EvtGen [12]. The generated particles are traced through the detector with 
GEANT4 [13], taking into account the details of the geometry and material composition of 
the detector. 

5° and B^ mesons are reconstructed from a selected meson combined with a 
vector particle (p° or K*'^). The selection requirements are kept as similar as possible for 
B^ —7- D^K*^ and B^ — t- D^p^. The four charged particles in the decay are each required 
to have a transverse momentum px > 300 MeV/ c for the daughters of the vector particle 
and pt > 250 MeV/ c (400 MeV/ c) for the pion (kaon) from the meson decay. The 

of the track impact parameter with respect to any primary vertex is required to be 
greater than 4. A cut on the absolute value of the cosine of the helicity angle of the vector 
meson greater than 0.4 is applied. The tracks of the meson daughters are combined to 
form a vertex with a goodness of fit x^/ndf smaller than 5. The B meson vertex formed 
by the and the tracks of the V meson daughters is required to satisfy x^/ndf < 4. 
The smallest impact parameter of the B meson with respect to all the primary vertices is 
required to be smaller than 9 and defines uniquely the primary vertex associated to the B 
meson. Since the B^ or B^ should point towards the primary vertex, the angle between 
the B momentum and the B line of flight defined by the line between the B vertex and 
the primary vertex is required to be less than 10 mrad. Finally, since the measured z 
position (along the beam direction) of the D vertex (zo) is not expected to be situated 
significantly upstream of the z position of the vector particle vertex (zy), a requirement 

of {zu — zy)! ^ a1 £1 + 0"^ y > —2 is applied, where a^, d and a^, y are the uncertainties 

on the z positions of the D and V vertices respectively. 

The selection criteria for the V candidates introduce some differences between the 
signal and normalisation channel due to the particle identification (PID) and mass win- 
dow requirements. The K*^ (p°) reconstructed mass is required to be within 50 MeV/c^ 
(150 MeV/c^) of its nominal value [14]. The selection criteria for the and vector mesons 
include identifying kaon and pion candidates using the RICH system. This analysis uses 
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the comparison between the kaon and pion hypotheses, DLLktt, which represents the 
difference in logarithms of hkehhoods for the K with respect to the tt hypothesis. The 
particle identification requirements for both kaon and pion hypotheses have been opti- 
mized on data. The thresholds are set at DLLktt > and DLLktt < 4, respectively, for 
the kaon and the pion from the . The misidentification rate is kept low by setting the 
thresholds for the vector meson daughters to DLLktt > 3 and DLLktt < 3 for the kaon 
and pion respectively. In order to remove the potential backgrounds due to 5^ — )■ Df7T~ 
and 5°— > D+TT" with D~ — K*^K~ and D~ — )■ K*°K~, vetoes around the nominal D~ 
and meson masses [14] of ±15 MeV/c^ are applied. Monte Carlo studies suggest that 
these vetoes are more than 99.5% efficient on the signal. 

Finally, multiple candidates in an event (about 5%) are removed by choosing the 
B candidate with the largest B flight distance significance and which lies in the mass 
windows of the and the vector meson resonance. 



3 Extraction of the ratio of branching fractions 

The ratio of branching fractions is calculated from the number of signal events in the two 
decay channels D°K*° and 5°-^ ^V", 

where the e parameters represent the total efficiencies, including acceptance, trigger, re- 
construction and selection, and fs/fd is the ratio of B^ and B^ hadronization fractions 
in pp collisions at y/s = 7 TeV. Since a given event can either be triggered by tracks 
from the signal or by tracks from the other B hadron decay, absolute efficiencies cannot 
be obtained with a great precision from the Monte Carlo simulation due to improper 
modelling of the generic B hadron decays. In order to reduce the systematic uncertainty 
related to the Monte Carlo simulation of the trigger, the data sample is divided into two 
categories: candidates that satisfy only the hadronic hardware trigger^ (TOSOnly, since 
they are Triggered On the Signal (TOS) exclusively and not on the rest of the event) and 
events which are Triggered by the rest of the event Independent of the Signal candidate B 
decay (TIS). Approximately 6% of candidates do not enter either of these two categories, 
and are vetoed in the analysis. The B^ — D^p^ signal yield is extracted separately for 
the two trigger categories TOSOnly and TIS; the B^ — )■ D^K*^ signal yield is extracted 
from the sum of both data samples. The ratio of efficiencies are sub-divided into the 
contributions arising from the selection requirements (including acceptance effects, but 
excluding PID), Vsei, the PID requirements, rpiD, and the trigger requirements, rTosoniy 

^Events passing only the muon trigger on the signal candidate tracks are rejected. 
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and Txis- The ratio of the branching fractions can therefore be expressed as 



-Slg. 



X — X Tsel X rpiD X 



B^-j. D°K*° 



(2) 



B (50 ^ DOpO) B {K*^ ^ K+TX-) f, 
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go ^ pOpO 
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where a represents a correction factor for the "non-p"" contribution in the — > D^p^ 
decays. 

The values of the efficiency ratios are measured using simulated events, except for 
rpiD = 1.09 ± 0.08 which is obtained from data using the D* — )■ D^n decay with 

—7- Ji~7r'^ where clean samples of kaons and pions can be obtained using a purely 
kinematic selection. Since the event selection is identical for the in the two channels 
of interest, many factors cancel out in r^ex = 0.784 ± 0.024 thereby reducing the system- 
atic uncertainties. The values of the trigger efficiency ratios, riosoniy = 1-20 ± 0.08 and 
Txis = 1-03 ±0.03, depend on the trigger configurations and are therefore computed from 
a luminosity-weighted average. The quoted uncertainties reflect the difference between 
data and Monte Carlo simulation mainly caused by the energy calibration of the trigger. 

The numbers of events in the two D^p^ trigger categories, -^^o^^^^o^o and ^^^pOpO, 
and A^^n' r-,n.^.n are extracted from a simultaneous unbinned maximum likelihood flt to 
the data. In order to simplify the description of the partially reconstructed background, 
the lower edge of the B meson mass window is restricted to 5.1 GeV/c^ for the i?*^— > D^p^ 
decay mode and to 5.19 GeV/c^ for the S° — >■ D^K*^ decay mode. There are four types 
of events in each category: signal, combinatorial background, partially reconstructed 
background and cross-feed.^ The signal B meson mass PDFs for B^ — )■ D^p^ and B^ — )■ 
D^K*^ are parametrized for each channel using the sum of two Gaussians sharing the same 
mean value. The mean and width of the core Gaussian describing the B^ — > mass 
distribution are allowed to vary in the flt. The fraction of events in the core Gaussian, 
0.81 ± 0.02, and the ratio of the tail and core Gaussian widths, 2.04 ± 0.05, are flxed 
to the values obtained from Monte Carlo simulation. In order to take into account the 
difference in mass resolution for the B^ — D^p^ and i?^ — D^K*^ decay modes, the value 
of the ratio of core Gaussian widths = 0.89 ± 0.03 is flxed from the Monte Carlo 



simulation. The mass difference between the means of the B^ and signals is flxed to 
the nominal value [14]. 

The combinatorial background mass distribution is modelled by a flat PDF and the 
partially reconstructed background is parametrized by an exponential function; the expo- 
nential slope is different in the B^ — > D^p^ and 5° — )■ D^K*^ categories. Since the number 
of S"— )■ D^p^ decays is larger than that of -8°— >■ D^K*^, the contribution from misiden- 
tifled pions as kaons from real 5°— > D^p^ has to be taken into account. The fractions of 
the cross-feed events, fDOpO^DOK*o = 0.062 ± 0.031 and fDOK*o^DOpO = 0.095 ± 0.047, are 

"^Thc cross-feed events are due to particle misidentification on one of the vector daughters; some I?°p° 
events can be selected as D'^/v *° and vice versa. 
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constrained using the results from a Monte Carlo study corrected by the PID misidenti- 
cation rates measured in data. The PDF for the cross-feed is empirically parametrised 
by a Crystal Ball function [15], whose width and other parameters are taken from a fit to 
simulated events in which — D^K*^ events are misidentified as — )■ D^p^ and vice 
versa; the width is fixed to 1.75 times the signal resolution. For the — > D^p^ decay 
mode, the events are further split according to the TOSOnly and TIS categories. 

In summary, 13 parameters are free in the fit. Four shape parameters are used, two for 
the signal and two for the partially reconstructed backgrounds. In addition, nine event 
yields are extracted, three (signal, combinatorial and partially reconstructed backgrounds) 
in each of the three categories: B°^ D^p° (TGSOnly and TIS) and D°K*°. 

The results of the fit for D^p^ and B^ — )■ D^K*^ are shown in Fig. 1 and Fig. 2. 
The overall signal yields are 154.1 ± 15.1 and 34.4 ± 6.8 respectively. The yields for the 
different components are summarised in Table 1. 




LHCb 

DV°, TOSOnly J 



5.3 5.4 5.5 5.6 



5.4 5.5 

m(DY) (GeV/c^) 




5.4 5.5 

m(D°p°) (GeV/c 



5.6 



Figure 1: The invariant mass distribution for the )■ D^p^ decay mode for the TOSOnly 
(left) and TIS (right) trigger categories with the result of the fit superimposed. The black 
points correspond to the data and the fit result is represented as a solid line. The signal is 
fitted with a double Gaussian (dashed line), the partially reconstructed background with 
an exponential function (light grey area) and the combinatorial background with a flat 
distribution (dark grey area) as explained in the text. The contributions from cross-feed 
are too small to be visible. 



Table 1: Summary of the fitted yields for the different categories. The background yields 
are quoted for the full mass regions. 



Decay mode 


Signal yield 


Part. rec. bkgd yield 


Comb, bkgd yield 




34.4 ±6.8 


17.5 ± 11.4 


29.8 ±8.4 


B^^ D^p^ (TOSOnly) 


77.0 ± 10.1 


55.4 ± 10.1 


95.5 ± 13.1 


(TIS) 


77.1 ± 11.2 


85.6 ±12.9 


176.0 ± 17.5 



In order to check the existence of other contributions under the vector mass peaks. 
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Figure 2: The invariant mass distribution for the — )■ D^K*^ decay mode with the 
result of the fit superimposed. The black points correspond to the data and the fit result 
is represented as a solid line. The signal is fitted with a double Gaussian (dashed line), 
the partially reconstructed background with an exponential function (light grey area), 
the combinatorial background with a flat distribution (dark grey area) and the cross-feed 
from B'^ — > D^p^ (intermediate grey area) as explained in the text. 



the sPlot technique [16] has been used to obtain background subtracted invariant mass 
distributions. The sWeights are calculated from the reconstructed B invariant mass dis- 
tribution using the same parametrization as in the analysis, the selection being the same 
except for the V invariant mass ranges which are widened. It was checked that there is 
no correlation between the B and the V invariant mass. The resulting plots are shown 
in Fig. 3, where the resonant component is fitted with a Breit-Wigner convoluted with 
a Gaussian and the non-resonant part with a second order polynomial. While the K*^ 
region shows no sign of an extra contribution, the p° region shows a more complicated 
structure. An effective "non-p°" contribution is estimated using a second-order polyno- 
mial: 30.1 ± 7.9 events contribute in the p*^ mass window (±150 MeV/c^). The measured 
7- D^p^ yields are corrected by a factor a = 0.805 ± 0.054 (see Eq. 2), consistent with 

expectations based on previous studies of the D^7t^it~ Dalitz plot [17, 18]. 

13( bo^dok*°) 

The ratio of branching fractions, \1 — / , is calculated using the measured yields 

of the B^ —7- D^p^ signal in the two trigger categories, corrected for the "non-p°" events and 
assumed to contribute proportionally to the TOSOnly and TIS samples, the 5^— !■ D^K*^ 

yield and the values of the r ratios quoted above. The result is ^(^^o^ijOpO^ ~ 1-48 ±0.34, 
where the uncertainty is statistical only. The small statistical correlation between the two 
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Figure 3: The p° (on the left) and K*^ (on the right) invariant mass distributions obtained 
from data using an sPlot technique. The level of non K*^ combinations in the — D^K*^ 
peak is negligible. Despite being mainly due to D^p^ combinations, the — t- D^p^ 
contains a significant contribution of "non-p°" events. The black points correspond to the 
data and the fit result is represented as a solid line. The resonant component is fitted 
with a Breit-Wigner convoluted with a Gaussian (dashed line) and the non- resonant part, 
if present, with a second-order polynomial (grey area). 



yields due to the cross-feed has been neglected. 



4 Systematic uncertainties 

A summary of the contributions to the systematic uncertainty is given in Table 2. The 
PID performances are determined with a D* — )■ D^n data calibration sample reweighted 
according to the kinematical properties of our signals obtained from Monte Carlo sim- 
ulation. The systematic uncertainty has been assigned using the kinematical distribu- 
tions directly obtained from the data. However, due to the small signal yield in the B^ 
case, this systematic uncertainty suffers from large statistical fluctuations which directly 
translate into a large systematic uncertainty on the kaon identification. The statistical 
uncertainty obtained on the number of "non-p°" events present in the p° the mass win- 
dow (±150 MeV/c^) has been propagated in the systematic uncertainty. The differences 
observed between Monte Carlo simulation and data on the values of the D'^ and vec- 
tor mesons reconstructed masses, as well as on the transverse momentum spectra, have 
been propagated into the uncertainty quoted on Vsci- The relative abundances of TOSOnly 
and TIS triggered events determined from simulated signal are in good agreement with 
those measured from data. This provides confidence in the description of the trigger in 
the Monte Carlo simulation. Since these relative abundances are directly measured in 
data, they do not enter the systematic uncertainty evaluation. However, the difference 
in trigger efficiency between the 5° — )> D'^p^ and the 5° D^K*^ decay modes is taken 
from Monte Carlo simulation; this is considered reliable since the difference arises due 
to the kinematical properties of the decays which are well modelled in the simulation. 



7 



The difference in the energy measurement between the hardware trigger clustering and 
the offfine reconstruction clustering is conservatively taken as a systematic uncertainty 
due to the hadronic trigger threshold. The systematic uncertainty due to the TIS trigger 
performances on the two decay modes is obtained assuming that it does not depend on 
the decay mode (rxis = !)• 

The systematic uncertainty due to the PDF parametrizations has been evaluated using 
toy Monte Carlo simulations where the different types of background have been gener- 
ated using an alternative parametrization (wide Gaussians for the partially reconstructed 
backgrounds, first order polynomial for the combinatorial backgrounds) but fitted with 
the default PDFs. 

The total systematic uncertainty is obtained by combining all sources in quadrature. 
The dominant sources of systematic uncertainty are of statistical nature and will be 
reduced with more data. The error on the ratio of the fragmentation fractions [10] is 
quoted ClS Sb separate systematic uncertainty. 



Table 2: Summary of the contributions to the systematic uncertainties. The uncertainty 


on the r ratio gives the range used for the systematic uncertainty extraction on the ratios 


of the branching fractions. 




Source 


Relative uncertainty 


Difference between data and MC to compute rpiD = 1.09 ± 0.06 


5.8% 


Uncertainty on the "non-p°" component a = 0.805 ± 0.054 


6.8% 


MC selection efficiencies r^^i = 0.784 ± 0.024 


3.1% 


LO Hadron threshold rxosoniy = 1.20 ± 0.08 


3.0% 


TIS triggering efficiency rjis = 1-03 ± 0.03 


1.6% 


PDF parametrisations 


1.0% 


Overall relative systematic uncertainty 


10.2 % 


Fragmentation fractions 


7.9% 



5 Summary 

A signal of 34.4 ± 6.8 -B°— D^K*^ events is observed for the first time. The significance 
of the background fluctuating to form the 5^ signal corresponds to approximately nine 
standard deviations, as determined from the change in twice the natural logarithm of 
the likelihood of the fit without signal. Although this significance includes the statistical 
uncertainty only, the result is unchanged if the small sources of systematic error that 
affect the yields are included. The branching fraction for this decay is measured relative 
to that for ^ D^p^, after correcting for the "non-p°" component, to be 

-^-^ = 1.48 ±0.34 ±0.15 ±0.12, 3 

13 {B^^ D^pO) ^ ' 
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where the first uncertainty is statistical, the second systematic and the third is due to the 
uncertainty in the hadronisation fraction (fs/fd)- 

The result is in agreement with other measurements of similar ratios and supports 
the SU(3) breaking observation in colour suppressed B^^s) ~^ decays. Using 

B(^B^^ D^P^^ = (3.2 ± 0.5) X 10~^ [14] for the branching fraction of the normalising 
decay, a measurement of the i?^ — D^K*^ branching fraction, 

B ^ = (4.72 ± 1.07 ± 0.48 ± 0.37 ± 0.74) x 10-^ (4) 

is obtained, where the first uncertainty is statistical, the second systematic, the third 
due to the uncertainty in the hadronisation fraction {fs/fd) and the last is due to the 
uncertainty of the — D^p^ branching fraction. A future, larger data sample will allow 
the use of the B^ — )■ D^K*^ decay as the normalising channel, which will reduce the 
systematic uncertainty. 
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